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1 Problem 

A typical high-energy neutrino beam is made from the decay of vr mesons that have been 
produced in proton interactions on a target, as sketched in the figure below. 
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^ I Suppose that only positively charged particles are collected by the "horn". The main 

m 
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^ ■ source of neutrinos is then the decay 7r+ —>■ ix'^v^. 



1. Give a simple estimate of the relative number of other types of neutrinos than in 
the beam. 

^ , 2. If the decay pions have energy E.,^ ^ ttItt, what is the characteristic angle 9c of the 

I decay neutrinos with respect to the direction of the tt"*"? 

^1 3. If a neutrino is produced with energy E^, ^ rriT^, what is the maximum angle 9niax{E,y) 

■ between it and the direction of its parent pion (which can have any energy)? What is 

the maximum energy Ei, at which a neutrino can be produced in the decay of a pion 
^ . if it appears at a given angle 6 with respect to the pion's direction? 

d • Parts 4 and 6 explore consequences of the existence of these maxima. 

4. Deduce an analytic expression for the energy-angle spectrum (PN / dEi^dVL for neutrinos 
produced at angle ^ < ^c* to the proton beam. You may suppose that E^, ^ m^, that 
the pions are produced with an energy spectrum dN/dE^^ oc {Ep — E^^Y , where Ep is 
the energy of the proton beam, and that the "horn" makes all pion momenta parallel 
to that of the proton beam. 

5. At what energy Ey^^cak does the neutrino spectrum peak for ^ = 0? 

6. Compare the characteristics of a neutrino beam at ^ = with an off- axis beam at angle 
9 such that -E'i/,max(^) is less than E^^^^akiP = 0). 

Facts: = 139.6 MeV/c^, r„ = 26 ns, = 105.7 MeV/c^, = 2.2 /is. In this 
problem, neutrinos can be taken as massless. 
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2 Solution 



In this solution we use units where c = 1. 

1. Besides the from the decay tt"*" — > ^Ji^v^-, the beam will also contain and z/g 
from the subsequent decay yU.^ e^z/gJ^^i- Both of these decays occur (primarily) in 
the "decay pipe" shown in the figure. As both the pions and muons of relevance are 
relativistic in this problem, they both have about the same amount of time to decay 
before they are absorbed in the "dump" . Hence, the ratio of number of muon decays 
to pion decays is roughly the same as the ration of there lifetimes, i.e., about 0.01. 
Our simple estimate is therefore. 



0.01. 



Experts may note that an additional source of Ve is the decay vr^ e^Ve at the level 
of 10~^. Also, mesons will be produced by the primary proton interaction at a rate 
about 10% that of vr"*". About 65% of decays are to n'^f^, which add to the main 
z/^ beam, but about 5% of the decays are to vr+vr^z/e, which increases the z/g component 
of the neutrino beam by about 0.1 x 0.05 = 0.005. 

2. Parts 2-6 of this problem are based on the kinematics of charged pion decay, which are 
closely related to kinematic features of neutral pion decay, 7r° — i- 77 |]1|. 

Experts may guess that the characteristic angle of the decay neutrinos with respect to 
the parent pion is = I/Ttt = ""^ttZ-E'tt- The details of the derivation are needed in 
part 3. 

We consider the decay ir ^ fiu in the rest frame of the pion (in which quantities will 
be labeled with the superscript and transform the results to the lab frame. 

Energy- momentum conservation can be written as the 4-vector relation, 

vr = /i + I/, (2) 

where the squares of the 4- vectors are the particle masses, vr^ = m^, /i^ = and 
z/^ = 0. As we are not concerned with details of the muon, it is convenient to rewrite 
eq. d) as 

/i = vr - z/, (3) 

and square this to find 

ml = ml- 2(71 -u). (4) 
In the rest frame of the pion, its 4-vector can be written 

7r = (m,, 0,0,0). (5) 

Taking the z axis to be the direction of the pion in the lab frame, the 4-vector of the 
(massless) neutrino in the pion rest frame can be written as 

z/ = (E*,E*sinr,0,E*cosr), (6) 
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since the energy and momentum of a massless particle are equal. The 4-vector product 
(tt • I/) = ttqUq — TTiVi is therefore 

(vr • i.) = m^El. (7) 
Hence, from eq. the energy of the neutrino in the pion rest frame is 

2 2 
ffi — TfL 

El = — ^ = 29.8 MeV, (8) 

using the stated facts. 

We can now transform the neutrino 4-vector (^) to the lab frame, using the Lorentz 

boost = -Evr/mTT, 

p = {Eu,EySm6,Q,EyCos6) 

= (7.^:(1 + cos r). El sin r, 0, ^.El{f3^ + cos 9*)). (9) 

The pion has spin zero, so the decay is isotropic in the pion rest frame. A relation for 
the angle 9 between the neutrino and its parent pion can be obtained from the 1 and 
3 components of eq. @, 

ElsiYi9* , , 

tan 9 = — — ^ — . 10 

The characteristic angle of the decay in the lab frame is usefully associated with decays 
at 9* = 90° in the pion rest frame. Thus, 

tan^c = ^- (11) 
When E.„ ^ niTr then 7^ ^ 1, ^ 1, and 

3. We now consider the lab angle ([l0|) between the neutrino and its parent pion with 
emphasis on the neutrino energy rather than the pion energy. If Ey ^ m^, then 
Et, ^ m^r also, so 7^ ^ 1 and ~ 1. Then we can write 

K*sin^* K>in^^ 

tan^?^ ^ , (13) 

7.^,^(1 + cos ^'^) E, ' ^ ^ 



using the time component of eq. (^. Since sin 9* cannot exceed unity, we see that there 
is a maximum lab angle 9 relative to the direction of the pion at which a neutrino of 
energy Ey can appear, namely 

^ El ^ 30 MeV 

C'max ~ „ ~ „ , \^^) 

which is small for <C E^. 

If instead, the angle 9 is given, eq. (0) also tells us that 

tant^ tant^ 
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4. We desire the neutrino spectrum in terms of the laboratory quantities E^, 6 and 0. 
We expect that the spectrum is uniform in the azimuthal angle 0. We are given the 
energy spectrum dN/dE^, oc {Ep — E^,)^ of the parent pions, and we have deduced that 
the spectrum is isotropic in the pion rest frame, i.e., flat in cos^*. Hence, we seek the 
transformation 



d^N d'^N 
oc 



d^N 



dE^dQ dEyd cos d dEj^d cos 9* 
where the Jacobian is given by 



.J (E^, cos 9* ;E^, cos 9) cx {E^-E^fj, (16) 



J{Et,, cos 9*] El,, cos( 



OEtt d cos 9* 
OEtt d cos 9* 



d cos 9 d cos t 



(17) 



The "exact" form of the Jacobian is somewhat lengthy, so we will simplify to the extent 
we can by noting that when E^, ^ myr, the parent pion has i?7r 3> rriT^ also, and so 
/^TT ~ 1- Also, part 3 tells us that 9 is very small for any value of 9*. 

We already have relation (|T3p between E,^, tan 6* and sin 6'*, so we can write 



cosr 



1 - sin^ 9* 



E^ 

% tan2 9 



\ 



El 



1 



E*2 VC0S2^ 



Thus, 



for small 9, and 



d cos 9* 
dE, 



d cos 9* 
dcos9 



J& tan^ 9 



|l tan^ 9 



_B;2 cosS 9 

l-^tan^9 



EJ' 



COS^* ' 



E?, 



Et? cos 9* 



We can also use time components of eq. (^ to write 



Hence, 



Et, 



E,. 



E,. 



E*{1 + (3^cos9*) E*{l + cos9*) 



(19) 



(20) 



:2i] 



dE,, 



and 



K*(l + cos^* 



m^Ey 



d cos 9* 



Et{l + cos9*)^ dE^ 
rrii^Ey d cos 9* 



E,, 



El9^ 



m^E* cos 9* 



ElE. 



dE^ 

9 COS 61 E*{\ + cos9*Y d cos9 m.„E*cos9* 



[22] 



(23) 
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The Jacobian (|T^) is therefore 



J 



+ 



m-TT E* cos ( 



E,e^ 

Ej2 cos e* 
El 



Et? cos 9* 



(24) 



m,r E* cos 0* E*'^cos6* 

and hence the neutrino spectrum can be written from eq. (0) as 



dEyd cos 6 



oc {Ep - E^ 



cos 9* 



(25) 



Because the factor cos 6* in the denominator of the Jacobian can go to zero, it is 
possible that the neutrino flux is higher for nonzero values of the lab angle 6. 

5. On the axis, 6 = 0, 6* = 0, and = mj,Ej,/2E* ^ 2E^ according to eq. (|2ll) . In this 
case, the neutrino spectrum (p5|) is 



d^Nje = 0) 
dEi^d cos 9 

The peak of the spectrum occurs at 



oc Er, 



2E* 



Et. 



E, 



^E 
7m^ ^ 8 



(26) 



(27) 



6. For an off- axis neutrino beam (at a nonzero value of angle 9) we must evaluate the 
spectrum (|25D using relations (|T^) and (|2T|). This is readily done numerically. For 
example, a plot of the pion energy E^^ needed to produce a neutrino of energy Ej^ at 
various angles 9 is shown below. 



> 




As expected from part 3, we see that for a given angle 6*, there is a maximum possible 
neutrino energy, and as the neutrino energy approaches this value, a large range of 
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pion energies contributes to a small range of neutrino energies. This will result in an 
enhancement of the neutrino spectrum. If we desire the enhancement at a particular 
neutrino energy, we should look for the neutrinos close to the angle 6'max given in 
eq. ([l^) , which is independent of the proton/pion energy. 

A numerical evaluation of the neutrino spectrum (^3p for several values of angle 9 with 
respect to the proton/pion beam is shown below. 



o 

■i 3 

3 
(D 

z 

2 
> 

CD 

OC 1 







e = 0.03 






E„ = 12 GeV 












^ e = 0.025 
















^ e = 0.02 








9 = 0.015 




/ e = 0.01 










^ ' "r- 




1 \ ' ' 



(GeV) 



We see that the spectrum of neutrinos at a nonzero angle is peaked at a lower energy, 
and is narrower, than that at zero degrees, due to the existence of a maximum possible 



neutrino energy (15) in decays at a given angle to the direction of the parent pion. This 
effect is especially prominent when E^rna.x{9) ^ (30 MeV)/9 is less than -E',y,peak(^ = 0), 
as then there is a substantial rate of higher energy pions all of which decay into a 
narrow band of neutrino energies at this angle. 

The spectral narrowing of an off-axis neutrino beam remains in more complete calcu- 
lation 1^, that include the nonzero transverse momenta of the pions before and after 
passing through the "horn" , although the spectrum will not have such hard edges, and 
the favored angle-energy combination is ^ ~ (50 MeV)/E^. 

In sum, the existence of a maximum energy for neutrinos that decay at a given angle to 
their parent pions implies that many different pion energies contribute to the this neu- 
trino energy, which enhances the neutrino spectrum at this angle-energy combination, 
6 ^ (30-50 MeV)/E^. 
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